#pragma once

void LU_decompose(double **a, double **L, double **U);  // 法2.2

void print_matrix(double** phi, double** eta, double **xx, int n, char* buf);

void LU_row(double **a);  // 讲义算法2.4  行优先存储，L和U直接存储在原来的a矩阵中

void Gauss_solver(double **a, double *x, double *b, int  n);  // 讲义算法2.4  行优先存储，L和U直接存储在原来的a矩阵中

void Gauss_partial_pivoting(double **a, double *x, double *b, int n);  // 讲义算法2.8 // 部分选主元高斯消去法 

void print_vector(double *x, int n, char* buf);

void mat_inverse(double **a, double **inverse, int n);
